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ABSTRACT. 

An  eigenspace  assignment  approach  to  the  design  of  parameter  insensitive  control  laws  for 
linear  multivariable  systems  is  presented.  The  control  design  scheme  utilizes  flexibility  in 
eigenvector  assignments  to  reduce  control  system  sensitivity  to  changes  in  system  parameters.  The 
methods  involve  use  of  the  singular  value  decomposition  to  provide  an  exact  description  of  allowable 
eigenvectors  in  terms  of  a  minimum  number  of  design  parameters.  In  a  design  example,  the 
methods  are  applied  to  the  problem  of  symmetric  flutter  suppression  in  an  aeroelastic  vehicle.  In  this 
example  the  flutter  mode  is  sensitive  to  changes  in  dynamic  pressure  and  eigenspace  methods  are 
used  to  enhance  the  performance  of  a  stabilizing  minimum  energy /  linear  quadratic  regulator 
controller  and  associated  observer.  Results  indicate  that  the  methods  provide  feedback  control  laws 
that  make  stability  of  the  nominal  closed  loop  systems  insensitive  to  changes  in  dynamic  pressure. 


tPartially  supported  by  the  grant  ONR-N00014-86-K-0693  and  by  NASA  Contract 
No.NASl-18107  while  the  author  was  in  residence  at  ICASE,  NASA  Langley  Research  Center, 
Hampton  VA  23665-5225. 


AN  APPLICATION  OF  EIGENSPACE  METHODS 
TO  SYMMETRIC  FLUTTER  SUPPRESSION 


1.  INTRODUCTION. 

A  common  problem  in  flight  control  applications  is  to  design  feedback  and  observer  gain 
matrices  so  that  the  system  remains  stable  and  satisifies  given  performance  requirements  over  a  range 
of  flight  conditions.  A  basic  approach  is  to  fix  a  flight  condition  as  a  design  point  and  to  synthesize 
gain  matrices  so  that  stability  and  performance  requirements  are  met  at  the  design  point.Nearby 
flight  conditions  yield  evalution  models.  If  performance  requirements  are  not  met  at  the  evaluation 
points,  one  returns  to  the  design  model  and  adjusts  the  existing  design  to  meet  more  stringent 
requirements.  In  such  cases  it  frequently  occurs  that  certain  variables  used  for  performance 
evaluation  are  more  sensitive  to  changes  in  the  design  point  than  others  and  a  reduction  in  sensitivity 
can  be  explicitly  considered  in  the  design  process.  In  this  paper  eigenspace  methods  are  presented 
for  the  design  of  parameter  insensitive  control  laws.  The  methods  are  applied  to  the  problem  of 
symmetric  flutter  suppression  in  an  aeroelasdc  vehicle. 

The  ability  to  shape  the  fundamental  modes  of  a  system  by  modifying  selected  eigenvalues  and 
eigenvectors  underlies  the  appeal  of  eigenspace  methods  for  control  system  design.  The  basic 
freedoms  and  limitations  of  eigenstructure  assigment  methods  by  full  state  feedback  or  output 
feedback  are  presented  in  the  papers  by  Moore  [12],  by  Srinathkumar  [14],  by  Andry,  Shapiro  and 
Chung  [2],  and  by  Kautsky,  Nichols  and  Van  Dooren  [7].  The  utility  of  eigenstructure  assignment 
methods  in  observer  designs  is  discussed  in  the  work  of  Kazerooni  and  Houpt  [8].  In  these  methods 
feedback  and  observer  gain  matrices  are  determined  to  yield  a  desired  eigenstructure.  Stability  and 
transmission  considerations  motivate  the  location  of  desired  eigenvalues.whereas  the  ability  to  shape 
the  system  response  [2,  12],  to  enhance  system  performance  [  5],  to  reduce  system  sensitivity  [  15], 
or  to  design  robust  control  laws  [7,  8]  motivates  the  selection  of  desired  eigenvectors. 

In  this  paper,  vector  space  methods  are  presented  for  eigenspace  assignment  which  take  into 
consideration  eigenvalue  sensitivity  to  plant  parameter  variations  and  performance  constraints.  The 
procedure  is  formulated  as  a  constrained  optimization  problem  in  Section  2  along  with  a  description 
of  eigenspace  assignment  methods  for  the  determination  of  feedback  and  observer  gain  matrices. 
Formulation  of  the  procedure  as  a  constrained  optimization  problem  allows  explicit  consideration  of 
costs  associatied  with  control  effort  in  the  design  procedure.  Eigenvalue  sensitivity  can  be  expressed 
explicitly  in  terms  of  closed  loop  eigenvectors  and,  consequently,  the  freedom  available  in  assigning 
eigenvectors  can  be  directly  related  to  sensitivity  reduction.  In  Section  2,  the  singular  value 
decomposition  |9]  is  used  to  provide  a  basis  for  the  attainable  eigenvectors  associated  with  a  desired 
eigenvalue.  Thus  an  explicit  parameterization/coordinatization  of  the  attainable  eigenvectors  is 
obtained.  An  alternate  algorithm  for  computing  a  basis  for  the  attainable  eigenvectors  was  presented 
in  the  work  of  Porter  and  D'Azzo  (  13]  and,  more  recently,  by  Kautsky,  Nichols  and  Van  Dooren 
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[7].  Some  eigenvector  assignment  procedures  require  the  control  designer  to  specify  a  desired 
eigenvector  and  the  nearest  attainable  eigenvector  is  then  computed  [2],  [5],  The  approach 
advocated,  here,  involves  a  precise  description  of  the  attainable  eigenvectors  and  results  in  the 
display  of  the  design  freedom  in  terms  of  a  minimum  number  of  independent  parameters. 

A  flutter  suppresssion/gust  load  alleviation  problem  is  formulated  in  Section  3  and  the 
results  of  an  extensive  design  example  are  presented  in  Section  4.  In  this  example,  the  flutter  mode 
is  sensitive  to  changes  in  dynmaic  pressure  and  eigenspace  methods  are  used  to  enhance  the 
performance  properties  of  a  LQR/LTR  designed  compensator.  Results  indicate  that  the  methods 
provide  feedback  control  laws  which  make  the  stability  of  the  nominal  closed  loop  system 
insensitive  to  changes  in  dynamic  pressure. 

2.  PARAMETER  INSENSITIVE  CONTROL  SYSTEM  DESIGN. 

In  this  section  a  procedure  based  upon  eigenstructure  assignment  methods  for  the  design  of 
parameter  insensitive  control  laws  is  presented.  A  review  of  eigenstructure  assignment  methods 
useful  in  the  design  of  full  state  feedback  control  laws  and  state  observers  is  included.  In  this  paper 
the  singular  value  decomposition  is  used  to  provide  a  coordinatization  of  the  allowable  eigenvectors, 
which  arise  in  full  state  feedback  designs.  These  coordinates  become  the  design  parameters  for 
sensitivity  reduction. 

Eigenstructure  assignment.  Consider  the  linear  system 

dx/dt  =  Ax  +  Bu  +  Tq  (1) 

where  A,  B,  and  F  denote  matrices  of  appropriate  dimensions  and  x,  u,  rj  denote  state,  control,  and 
disturbance  variables  respectively.  Assuming  full  state  feedback  control,  i.  e.  u  =  -Kx  ,  the 
eigenvalues  of  A  -  BK  can  be  arbitrarily  assigned  through  a  proper  choice  of  the  gain  matrix  K  if 
and  only  if  the  system  (1)  is  completely  controllable.  If  (1)  is  uncontrollable  then  uncontrollable 
eigenvalues  cannot  be  altered  by  a  state  variable  feedback  control  law  u  =  -Kx. 

The  ability  to  assign  eigenvalues  and  eigenvectors  through  state  variable  feedback  control  laws 
is  summarized  as  follows  [3,  7, 12,  14].  Let  A  =  {  Xj :  i=l, ...,  n  }  be  a  self  conjugate  set  of 
distinct  complex  numbers.  Here  A  denotes  a  set  of  desired  closed  loop  eigenvalues  and  must  include 
all  uncontrollable  eigenvalues  of  the  system.  There  exists  a  feedback  gain  matrix  K  such  that  (A  - 
BK)  Vj  =  Xj  vj  for  i  =  1, ...,  n  if  and  only  if 
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v*v  f  c  1  if  in  >i  •  a  v  >.vtor  whose  components  are  the  coordinates  o<  I  v  .  kv )'  f  v.  ith  ro-quv'  '  A 
*vUi»ro>:  I  V. .  I  *;ii.>  i1')  provides  a  parameterization  or  co< »rc1i ni*t i ytat io< '  «'♦*  the  al’owahV 
eit:cn\cctors.  If  K  is  a  conunllahle  eigenvalue  it  may  he  shown  that  die  d  m-m  •  ol 
Ker  |  a.  1  A.  B  i  equals  the  number  of  columns  of  B  and  the  columns  ol  W  j .  sec  ( 3  i  abo\ e,  are 
independent  if  B  is  of  full  rank  1 12).  The  singular  value  decomposition  1 7.  H).  SVD,  provides  an 
efficient  numerical  method  to  compute  an  orthonormal  basis  for  Ker  |  X  1  -  A.  B  j  .  If  the  SVD  of  j 
/.  I  -  A.  B  |  is 

I  0 

|  A  I  -  A.  B  )  =  [  Up  U,  |  ()  ()  |  VrV2  |'  (4) 

then  the  columns  of  W  =  V2  form  an  orthonomial  basis  for  the  Ker  |  X  I  -  A,  B  | 

Equation  (3)  provides  a  parameterization  or  coordinatization  of  the  attainable  eigenvectors 
associated  with  X  and,  consequently,  completely  describes  the  freedom  in  assigning  a  particular 
eigenvector.  One  should  note  that  the  number  of  free  parameters  describing  an  attainable 
eigenvector  equals  the  rank(  Wj)  <  dim  Ker  |  X  I  -  A,  B  ] .  These  free  parameters  will  be  utilized 
as  control  design  parameters  in  the  application  to  follow.  Further  characterizations  of  the  attainable 
eigenvectors  and  the  number  of  degrees  of  freedom  in  eigenvector  assignment  are  presented  in  the 
work  of  Kautsky,  Nichols  and  Van  Dooren  (7). 

In  this  design  example  the  coordinates  of  an  attainable  eigenvector  are  modified  to  obtain 
feedback  control  laws  which  are  insensitive  to  changes  in  a  model  parameter.  The  approach  requires 
that  derivatives  of  eigenvalues  with  respect  to  system  parameters  be  computed.  It  is  well  known  |7. 

1 0 1  that  if  X  is  a  distinct  eigenvalue  of  a  matrix  Q(q),  depending  on  a  parameter  q,  and  if,  v  are  left 
and  right  eigenvectors  corresponding  to  X  ,  respectively,  then 

dX  u'dQ/dqv 
cl  c]  u'  v 

The  dependence  of  this  sensitivity  on  the  eigenvectors  u’  and  v  indicates  a  potential  reduction  in 
sensitivity  by  a  suitable  choice  of  eigenvectors.  Bounds  on  eigenvalue  sensitivities  are  directly 
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related  to  the  condition  number  of  the  modal  matrix  for  the  closed  loop  systems  and  are  discussed 
in  [7]. 

Eigenstructure  assignment  in  observer  designs.  In  a  recent  paper  [8],  Kazerooni  and  Houpt 
presented  a  procedure  for  loop  transfer  recovery  based  upon  eigenstructure  assignment  of  observers. 
This  procedure  for  observer  design  will  be  used  in  the  example  to  follow.  These  same  methods  have 
been  previously  applied  in  the  design  of  an  active  flutter  suppression  system  by  Garrard,  Liebst,  and 
Farm  [  5  ].  For  completeness  the  procedure  is  reviewed  in  this  section. 

Suppose  now  only  output  feedback  is  available  and  an  observer  is  to  be  employed  to  estimate 
the  state,  then  equation  (1)  is  coupled  with  the  following  output,  observer,  and  feedback  equations: 

y  =  Cx  (6) 

dz/dt  =  Az  +  Bu  +  H(y-Cz)  (7) 

u  =  -K  z  (8) 

Here  K  and  H  denote  feedback  and  observer  gain  matrices,  respectively. 

Setting  u  =  -Kz+  uQ  in  place  of  (8)  and  e  =  x  -  z,  one  obtains 

de/dt  =  (  A  -  HC)  e  +  B  u0  +  T  Tj .  (9) 

In  terms  of  transfer  functions  Tzu(s)  =  Txu(s)  -  Tgu(s).  For  recovery  it  is  desireable  to  have 
Tzu(s)  =Txu<s)  or  Teu(s>  =0- 

The  condition  Tgu(s)  =  0  can  be  described  in  terms  of  transmission  properties  of  (1)  and  (7). 
A  complex  number  £,  is  called  an  invariant  zero  [1 1]  of  the  system  (1)  and  (6)  with  left  zero  direction 
[v' ,  p’  1  provided 


[  V,  p*  ] 


^  I  -  A  B 
C  0 


and  [v' ,  p’  ]  *  0.  If  a  gain  matrix  H  can  be  chosen  which  assigns  the  eigenvalues  { Xj :  i=l,  ...,n) 
of  A  -  HC  to  invariant  zeros  {  ^  j :  i=l, ...,  n}  of  (1)  and  (6)  and  left  eigenvector  vj'  to 
corresponding  left  zero  direction  Vj'  then,  necessarily, 


I  -  A  B 


[v.,v.H] 


Moreover,  since  the  zero  state  response  of  (9)  is 

A- 1 

e(t )  =  2*  uj  vj'  B  e  J  *  uQ(t ) 
j=i 

where  uj  vj  denote  right  and  left  eigenvectors  of  A  corresponding  to  Xj,  it  follows  that  e(t)  =  0,  i.e. 
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Teu  =  0.  Thus  for  recovery  it  is  desirable  to  place  the  eigenvalues  of  A  -  HC  at  the  invariant  zeros 
of  (1)  and  (6)  and  the  left  eigenvectors  at  corresponding  left  zero  directions. 

Numerical  solution  of  (10)  reduces  to  the  solution  of  a  generalized  eigenvalue  problem  and 
existing  computer  software  [5, 6, 7]  can  be  used  to  obtain  the  invariant  zeros  and  left  zero  directions. 
In  this  procedure  the  desired  eigenstructure  represents  an  ideal.  There  is  no  guarantee  that  the 
desired  eigenstructure  is  attainable.  In  the  design  example  to  follow  nearest  attainable  eigenvectors 
are  computed  using  the  methods  presented  in  [6].  Once  desired  eigenvalues  and  attainable 
eigenvectors  have  been  specified,  the  observer  gain  matrix  can  be  computed  using  eigenstructure 
assignment  methods.  For  further  discussion  on  the  application  of  this  procedure  see  [  3]  and  [6].  If 
the  observer  gain  H(p)  depends  on  a  parameter  p  and  H(p)/p  — >  BW  as  p  —>  °°  where  W  is  a 
nonsingular  matrix  then  a  standard  arguement  [10]  implies  that  the  eigenvalues  of  A  -  HC  converge 
to  the  transmission  zeros  of  (1)  and  (6).  The  control  transfer  function  for  the  system  described  by 
equations  (1),  (6),  (7),  and  (8)  is  u  =  -  K(s)y  with 

K(s)  =  K  (  si  -  A  +  BK  +  HC  r1  H. 

Let  G(s)  denote  the  transfer  function  associated  with  (1)  and  (7),  i.e. 

G(s)  =  C(sl- 

In  this  case  Doyle  and  Stein  [6]  have  also  shown  that  K(s)G(s)  -»  KG(s)  pointwise.  Thus,  from 
another  viewpoint,  it  is  desireable  to  place  the  eigenvalues  of  the  observer  at  the  transmission  zeros 
of  the  system  (1)  and  (6). 

Parameter  insensitive  control  design.  Here  the  objective  is  to  determine  the  design  parameters  which 
minimize  the  sensitivity  of  certain  performance  variables  while,  at  the  same  time,  maintaining  other 
performance  variables  within  prescribed  bounds.  By  introducing  an  appropriate  penalty  function, 
this  problem  can  be  formulated  as  a  constrained  optimization  problem. 

In  this  paper,  the  design  variables  are  those  eigenvalues  and  eigenvectors  of  the  system 
matricies  A(qo)  -  B(  q<)  )K  and  A(  qc  )  -  HC(  q0  )  which  may  be  modified  in  order  to  reduce 
sensitivity.  The  design  parameters  are  the  real  and  imaginary  parts  of  the  designated  eigenvalues 
and  the  coordinates  of  the  associated  eigenvectors.  In  the  design  process  employed  herein,  a  basis 
for  the  allowable  eigenvectors  must  be  computed  for  each  eigenvector  to  be  modified,  see  equation 
(3).  The  coordinates  of  an  eigenvector  with  respect  to  this  basis  become  the  design  parameters  that 
relate  to  eigenvector  selection.  If  a  desired  eigenvalue  is  to  be  left  unaltered  throughout  a  design,  that 
is  eigensystem  design  freedoms  are  not  fully  utilized,  then  this  basis  need  only  be  computed  once. 

Let  a  denote  the  vector  of  design  parameters.  In  the  design  procedure  certain  performance 
variables,  those  which  measure  sensitivity  to  plant  variations,  are  to  be  minimized.  Let  s(  a)  denote 
a  vector  of  such  variables.  In  the  design  example  presented  in  Section  4,  s(  a)  will  denote 
eigenvalue  sensitivity  to  plant  parameter  variation.  Let  p(  a)  denote  a  vector  of  performance 
variables  to  be  kept  within  prescribed  bounds.  The  prescribed  lower  and  upper  bounds  for  the  i 
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Figure  3.3  Locus  of  Open  Loop  Poles  with  Dynamic  Pressure. 


At  the  design  point  the  controlled  system 

dx/dt  =  (A(  qj  -  B(  q0)K(  qc))  x  +  T(  eg  T] 

exhibits  good  performance  and  disturbance  rejection  properties.  The  rms  values  listed  in  column 
LQR  of  Table  4.1  indicate  acceptable  control  surface  activity  and  loads  for  a  12  ft./sec.  rms  gust 
velocity.  In  addition,  the  minimum  singular  value  for  the  return  difference  matrix  I  +  K(  qj 
(si- A(q0))-  ^ 3(^1  over  the  frequency  range  10  £  to  <,  300  is  approximately  one,  indicating  good 
disturbance  rejection  properties. 

At  this  design  point  an  observer  was  determined  by  attempting  to  place  the  eigenvalues  of  A 
-  HC  at  the  transmission  zeros  of  the  plant,  C(  qc)(  si  -  A(  q0))'^B(  qo),  and  corresponding 
eigenvectors  at  the  left  zero  directions.  The  open  loop  eigenvalues,  invariant  zeros  and  actual 
assignments  are  listed  in  Table  4.2.  Here  eigenvalues  1-10  correspond  to  position  and  rate 
variables,  11-15  correspond  to  unsteady  aerodynamic  states,  16-24  correspond  to  actuator  states, 
and  25-26  correspond  to  gust  states.  The  zeros  fall  in  three  categories:  those  with  large 
magnitude ,  those  of  the  same  order  of  magnitude  and  those  with  small  magnitude  relative  to  the 
magnitude  of  the  open  loop  eigenvalues.  The  six  small  magnitude  zeros  result  from  computattional 
errors  and  should  be  set  to  zero.  Actual  pole-zero  assignments  are  listed  in  columns  I  and  II  of 
Table  4.2.  The  desired  eigenvalues  and  eigenvectors  for  the  assignment  listed  in  column  I  of  Table 
4.2  were  made  as  follows:  zeros  of  the  same  order  of  magnitude  as  those  of  the  open  loop  poles 
(and  corresponding  left  zero  direction)  were  assigned  to  the  nearest  open  loop  eigenvalue  (and 


Control 

Deflections  (deg) 
Control 

Rates  (deg/sec) 


R/LTR 


Bending  Moment  (in  lbs)  26,344 

Shear  (lbs)  492.0 

Torque  (in  lbs)  1657. 


26,396 

492.9 

1708. 


Table  4.1  RMS  Performance  at  Design  Point  for  LQR  and  LQR/LTR  Compensators. 


Eigenvalue 

Zero 

I 

II 

1 

-13.8+430j 

-6.017+346.9j 

-6.017+346.9j 

-6.017+346.9j 

z 

3 

-19. 2+27 5j 

.1075 

-.12 

-4.8925 

4 

— 

-.1075 

-.1 

-5.1075 

5 

-1.4+194j 

-1.040+198.2j 

-1.040+198.2j 

-1. 040+198. 2j 

o 

7 

-43.6+ 136j 

.03479 

-.04 

-4.9652 

8 

— 

-.03479 

-.02 

-5.0348 

5.2+1 19.7j 

.000003+.  17j 

-.000003+.  17j 

-5.000003+.17j 

-92 

-111.2 

-111.2 

-111.2 

12 

-113.6 

-129.8 

-129.8 

-129.8 

13 

-112 

-122.2+.06059j 

-122.2+.06059j 

-122.2+.06059j 

14 

-112.4 

-122.2-.06050j 

-122.2-.06059j 

-122.2-.06059j 

15 

-112.2 

-111.7 

-111.7 

-111.7 

16 

-179.9 

oo 

-879.9 

-879.9 

17 

1  fi 

-125.6+287.9j 

oo 

-825.6+287.9j 

-825.6+287.9j 

lo 

19 

-175 

oo 

oo 

-875 

-875 

m 

-1 27.4+29 1.9j 

oo 

-827. 4+29 1.9j 

-827.4+29 1.9j 

21 

— 

oo 

— 

- 

22 

-185 

oo 

-885 

-885 

23 

-123.9+283.9j 

oo 

-823.9+283.9j 

-823.9+283.9J 

Z4 

25 

-.492 

oo 

-.492 

-.492 

-492 

26 

-.497 

-.497 

-.497 

-497 

Table  4.2  Eigenvalue-Zero  Assignments  for  Observer  Design. 
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v.  unstable.  Additional  calculations  show  that  the  matrices  A(q)  -  H(  q  )C  possess 
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ci;rc  4.2  1  (X  us  of  Pole  Locations  for  LQR  Compensators  A(q)  -  B(q)K(q  ) . 
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unstable  eigenvalues  for  q=  5.141  and  6.639  lb/in2-  It  is  desired  to  reduce  this  sensitivity  while 
maintaining  performance  requirements  by  use  of  the  eigenspace  assignment  methods  discussed  in 

previous  section. 

I  o  illustrate  the  design  procedure,  a  design  involving  only  a  small  number  of  parameters  is 
undertaken.  The  minimum  energy  stabilizing  controller  K(  qo)  modifies  the  eigenstructrure  of 
\f  q  i  in  a  simple  manner,  here  the  unstable  eigenvalues  are  flipped  symmetrically  with  respect  to 
: ne  imaginary'  axis  while  all  other  eigenvalue/eigenvector  pairs  remain  unchanged.  That  is  the 
stabilizing,  minimum  energy  feedback  control  law  only  alters  the  mode  shape  associated  with  the 
•pen  loop  instability.  In  the  design  example,  the  closed  loop  eigenvalues  are  left  in  the  stable 
m  vation  achieved  by  the  minimum  energy  controller  K(  qo)  and  an  attempt  is  made  to  reduce  the 
sensitivity  of  this  control  law  to  changes  in  dynamic  pressure  by  modifying  the  eigenvector 
. .  u  re-  ponding  to  the  unstable  mode.  Letting  X  denote  the  desired  eigenvalue  associated  with  the 
uist,:ble  pole,  the  attainable  eigenvectors  v  and  modified  gain  matrix  Kj^j  must  satisfy 


XI  A,  B  1 


■  s.mrpie  the  dimension  of  Ker[X  I-A,  B]  is  three.  Since  this  is  a  complex  vector  space  there 
•ix  tree  parameters  which  describe  the  attainable  eigenvectors  and  modified  gain  matrix 
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feedback  control  law  which  reduces  sensitivity  to  changes  in  dynamic  pressure,  maintains 
performance  requirements  and  is  robust  with  respect  to  model  uncertainties.  Although  all  of  these 
design  objectives  will  not  be  met,  the  example  does  illustrate  that  many  of  them  can  be  achieved  by 
simply  modifying  the  mode  shape  associated  with  the  unstable  eigenvalues. 

The  sensitivity  of  the  unstable  eigenvalue  to  changes  in  dynamic  pressure  is  calculated  from 
the  formula 

dX  /dq  =  {u-  (d(  A(q)  -  EKqjKCqJ)  /dq)  v  }  /  (u’  v) 

where  u'  and  v  denote  left  and  right  eigenvectors  of  A(q)  -  B(q)K(q0)  corresponding  to  the 
unstable  eigenvalue  and  d(  A(q)  -  B(q)K(qc))  /dq,  is  dX(  qQ)  /dq  is  estimated  by  the  difference 
quotient  { A(q)  -  B(q)K(qo)  -  (A(qo)  -  B(  q0  )K(qo)}/(q  -  qo).  The  sensitivity  is  dA.  /dq  =  43.88  + 
18.67j.  The  design  procedure  is  to  utilize  the  freedom  in  eigenvector  assignment  in  order  to 
minimize  the  magnitude  of  dX(  qc)  /dq  subject  to  the  stability  and  performance  constraints  at  this 
value  of  dynamic  pressure.  Robustness  constraints  are  not  explicitly  taken  into  consideration. 

Let  W  =  [  Wj',  W2'  ]'  be  a  basis  for  the  ker[  A.I  -  A,  B  ]  where  X  =  -5.2  +  1 19.7j.  The 
unstable  open  loop  eigenvalue  is  5.2  ±  1 19.7j.  The  attainable  eigenvectors  are  linear  combinations 
of  the  columns  of  W j ,  i.e.  v  =  Wj  c  for  an  arbitrary  vector  c.  The  design  parameters  a  =  c  are 
the  components  of  an  attainable  eigenvector  with  respect  to  the  basis  Wj,  which  is  obtained  using 
the  SVD  as  described  in  Section  2.  Here  s(  a)  =  d  X/  dq.  The  components  of  the  performance 
vector  p(  a)  are  a  stability  indicator,  position  and  rate  rms  values  for  each  controller,  and  wing 
bending  moment,  shear,  and  torque  rms  values.  The  stability  indicator  was  taken  taken  as  the 
maximum  of  the  real  part  of  the  closed  loop  system  eigenvalues.  In  our  design  example  the 
prescribed  lower  bounds  were  LB  =  (-10,000,0,0,0,0,0,0,0,0,0)  and  the  prescribed  upper  bounds 
were  UB  =  (0,3.7,15,15,372,372,372,30000,1000,2000).  The  weight  matrices  were  Qj  =  1  and 
Q2  =  diag(  10000,50,50,50,25,50,50,50,50,50). 

A  search  over  the  six  dimensional  parameter  space  determines  a  desired  eigenvector  to 
achieve  this  minimum  and  the  corresponding  gain  matrix  K^.  The  sensitivity  achieved  at  this  point 
is  dX(  q0)  /dq  =  10.16  +  2.07j.  The  corresponding  gain  matrix  is  listed  in  Table  A2  of  Appendix  A. 
To  further  illustrate  the  sensitivity  reduction  of  this  design  the  closed  loop  poles  for  the  systems 
A(q)  -  B(q)Kjvj  for  various  values  of  dynamic  pressure  have  been  depicted  in  Figure  4.3.  Note  that 
stability  is  achieved  at  each  evaluation  point.  Corresponding  rms  values  are  listed  in  Table  4.3.  It  is 
most  intersesting  that  these  results  can  be  achieved  by  simply  modifying  the  mode  shape 
corresponding  to  the  unstable  mode.  The  maximum  and  minimum  singular  values  of  the  return 
difference  matrix  I  +  Kj^  C(  q0  )  (si  -  A(  qc)  j'1  B(  qQ)  over  a  specified  frequency  range  have 
been  graphed  in  Figure  4.4.  This  figure  indicates  that  the  robustness  properties  of  the  minimum 
energy  controller  have  not  been  maintained  and  that  a  more  reasonable  design  would  include 
singular  value  constraints. 


Figure  4.3  Locus  of  Pole  Locations  for  Modified  LQR  Compensators  A(q)  *  B(q)Kj^j 


60  30 


rr.'*  | 


I 


10  100  Frequency 

Figure  4.4.  Maximum  and  Minimum  Singular  Values  of  the  Return 
Difference  Matrix  I  +  KM(sI  -  A)  1 B. 


Although  the  separation  principle  guarantees  the  stability  of  the  compui>atiT  .•  itS;  u \u  r. 
gain  matrix  K  =  Kyj  and  observer  gain  matrix  H  -  1 1(  q  ) .  instabilities  still  arise  at  off  design 
points  due  to  use  of  the  observer  gain  matrix  H  =  H(  q  ).  Consequently  a  modified  observer  dc.Mg;. 
was  attempted,  the  desired  pole  locations  for  this  modi  bed  design  are  listed  in  column  11  of  1  .-.bh 
4.2  and  the  resulting  gain  matrix  is  listed  in  Table  A4  of  Appendix  A  .  The  overall 
performance  of  the  compensator  may  now  be  evaluated.  First,  it  should  be  noted  that  numerical 
computations  verify  the  stability  of  the  compensator  at  all  evaluation  points.  Thus  the  design 
procedure,  i.e.  modification  of  the  eigenvectors  associated  with  the  unstable  mode  at  the  design 
point  for  feedback  gain  matrix  design  and  assignment  of  observer  eigenstru  rturc  to  corresponding 
system  zeros  and  invariant  directions  for  observer  gain  matrix  design  has  resulted  in  a 
compensator  whose  stability  is  less  sensitive  to  changes  in  dynamic  pressure.  However  one  should 
also  note  that  the  mis  values  for  the  system  with  compensator,  see  Table  4.4,  indicate  a  high  level 
of  control  activity  at  the  evaluation  points  furthest  from  the  design  point.  Further  improvement  in 
performance  should  be  attainable  though  modification  of  other  mode  shapes  and  by  adjusting  the 
weight  matrices  in  the  design  procedure  to  obtain  an  acceptable  blend  between  sensitivity  reduction 
and  control  activity.  Plots  of  the  maximum  and  minimum  singular  values  of  the  return  difference 
matrix  I  +  (si  -  A  +  BK^  +  H^C)" 1  Hj^C(sI  -  A)*^B  are  depicted  in  Figure  4.5.  Again, 
robustness  properties  of  the  minimum  energy  controller  have  not  been  retained. 


Dynamic  Pressure 


RMS  ^ 

4.4 17f  lh/in2 1 

4  476K 

5.141 

5.5.37 

6.639 

C  nnfrol 

3.008 

3.716 

4.989 

5.587 

6.506 

Deflections 

.380 

.4394 

.5725 

.6444 

.9127 

(deg) 

2.431 

3.026 

4.061 

4.518 

4.850 

C  nnrrol 

376.0 

465.4 

609.8 

659.9 

350.8 

Rates 

42.66 

52.0 

67.10 

71.38 

35.83 

(deg/sec) 

307.4 

380.6 

^  498.4 

538.6 

285.9 

Bending 

Moment  (in  lb) 
Shear  (lbs) 

24,825 

453.8 

458.2 

22,250 

484.3 

616.6 

27,950 

521.9 

848.5 

29,910 

563.7 

970.4 

39,008 

745.3 

2143 

Torque  (in  lbs) 

i 
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Table  4.3  RMS  Performance  for  Modified  LQR  Compensator  dx/dt  =  ( A(q)-B(q)K^)x. 
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F  igure  4.5.  Maximum  and  Minimum  Singular  Values  of  the  Return 

Difference  Matrix  1  +  KM  (si  -  A  +  BKM  +  HM  C)1  HmC(sI-A)  'b. 
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CONCLUDING  REMARKS. 

A  procedure  based  upon  eigenstructure  assignment  to  reduce  control  system  sensitivity  to 
changes  in  system  parameters  has  been  introduced.  Through  an  extensive  design  example  it  has 
been  shown  that  modification  of  a  single  mode  shape  can  lead  to  a  gain  matrix  whose  sensitivity  to 
parameter  variation  is  significantly  reduced.  The  design  example  also  indicates  that  a  proper  blend 
must  be  achieved  between  sensitivity  reduction  and  performance  constraints.  Another  problem  that 
needs  to  be  addressed  is  the  explicit  incorporation  of  robustness  requirements  into  the  design 
procedure.  Current  effort  involves  the  incorporation  of  explicit  robustness  constraints  in  the 
design.  Additional  freedoms  in  eigenvector  assignment  are  being  used  to  achieve  this  objective. 
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Appendix  A 
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